Macroscopic Chirality Fluctuations in Heavy Ion Collisions 
should induce CP forbidden Decays 
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If large fluctuations of quark chirality occur in heavy ion collisions, they result in macroscopic 
CP-odd "spots" of the so called theta- vacua, with a non-zero 8(x). We consider particular decays 
of mesons, CP-forbidden in the vacuum with zero 9, like 77 — > nn. We evaluate their rates for such 
decays near hadronic freezeout. These rates, as well as charge asymmetries already observed, are 
proportional to square of the CP-violating parameter (9 2 ) averaged over the fireball and events. 
With such input, we found that the forbidden decay rates are likely to be orders of magnitude larger 
than CP-allowed ones. We further estimated that up to about one per mill of r\ mesons produced 
in heavy ion collisions should decay in this way. We further discuss how those can be observed. 
We argue using STAR data on charge asymmetries for AuAu and CuCu collisions that the size of 
CP-odd spots at freezeout is as large as Cu nuclei: this fortunate fact (not explained so far by itself) 
suggests that the two-pion invariant mass is modified by only about a percent, which is comparable 
to experimental resolution. If so, we think experimental observation of these decays is within the 
reach of current dataset. If those decays are found, it would confirm that CP-odd interpretation of 
charge asyymetry is correct, even without complication related to geometry, impact parameter or 
magnetic field induced on the fireball. 
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I. INTRODUCTION 

Discovery of instantons [1 have revealed the dynami- 
cal role of gauge field topology, in particular existence of 
the #-vacua and the so called strong CP problem. With- 
out going into discussion of its cosmological aspects and 
axions, let us only mention how it plays out in a "Little 
Bang" , the heavy ion collisions. One notable suggestion, 
by Kharzeev, Pisarsky and Tytgat [5], was that analogs 
of theta- vacua can possibly appear at finite temperature, 
as metastable CP-odd bubbles. This might happened if 
the rj' mass can get very small near T c . However, lat- 
tice measurements of the topological susceptibility near 
T c , such as e.g. [TO], seem to exclude this scenario. But 
there is no need to have metastability of the bubbles: in 
fact initial state fluctuations of chirality would evolve in 
diffusively, as it is preserved in the QGP. The strongly 
coupled nature of it, manifesting itself in small viscosity 
and diffusion constant for charm, only helps to keep the 
CP-odd spot localized. 

But before we get to topological effects in heavy ion 
collisions, let us discuss the QCD vacuum first. The first 
empirical information on instantons in the QCD vacuum 
has been provided long ago by one of us j2] : based on the 
values of gluon and quark condensates it was conjectured 
that the mean size and the density (instantons plus anti- 
instantons together) are 

Pmst ~ 0.3 /m, n inst f» 1 fmr 4 (1) 

These numbers were confirmed by direct lattice obser- 



vations a decade later. For details of that and for sum- 
mary on QCD instanton-induced effects in chiral symme- 
try breaking and hadronic spectroscopy see [3]. 

The role of instantons in high energy collisions have 
been discussed in [4], [5], in which it was suggested that in- 
stantons are responsible for the so called "soft pomeron" , 
operating at a momentum scale ~ 1/p ~ IGeV, less 
than that for perturbative BFKL Pomeron. Eventually 
it lead to explanation of how topology/chirality fluctua- 
tions should occur in high energy collisions. It has been 
found that instantons, being perturbed by a collision, 
are excited and explode into certain states of glue and 
fermions which can be described semiclassically. Transi- 
tion from virtual (Euclidean) to real (Minkowskian) part 
of the semiclassical path proceed via the so called "turn- 
ing states" , at which the momentum (electric field) van- 
ishes . Such "turning states" have been identified as the 
so called Carter-Ostrovsky-Shuryak (COS) sphalcrons 
|6J. Analitic solutions for them, as well as classical de- 
scription of their subsequent decay have been made in the 
same paper. Further clarification of the fermion produc- 
tion in sphaleron explosion has been made by Shuryak 
and Zahed in [7 ], by solving Dirac equation in the ex- 
ploding sphaleron background. It clearly shows how chi- 
rally asymmetric set of quarks unULdudLSRSL are pulled 
from the Dirac sea into the positive energy states, violat- 
ing chirality by 6 units (as anomaly/index theorem equa- 
tion predicts) , with the calculated spectrum. It has been 
proposed that one can identify and experimentally study 
COS sphaleron production in double-diffractive pp col- 
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lisions, and currently there is some proposal to perform 
such experiments at RHIC. 

Returning to Au-Au collisions, one may use these ideas 
in order to estimate the number of sphalerons produced 
[8]. If we only consider the number of vacuum instantons 
in an initial state "pancake" then we get a so minimal 
number 



N sphalerons ^ ^instPinst^-^. 
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The maximal number in the whole fireball is perhaps an 
order of magnitude larger [8] . Assuming random statis- 
tics, this implies that the fluctuations of the topological 
charge in the the fireball is AQ top ~ y/ N spha i erons = 
3 — 10. As the size of the sphalerons is ~ p <C Rau, 
each of them can be viewd as exploding shells. In QGP 
those evolve diffusively and merge, in a process not yet 
understood quantitatively, producing at the time of fire- 
ball freezeout some chirally asymmetric spots, such as 
one shown in Fig.Q. 

Further development was related with the presence of 
a magnetic field B: the so called Witten effect states 
that when the nonzero 9 is present, it mixes with the 
electric field E ~ 9B. This leads to a current and even- 
tually charge separation in the direction of B . The mi- 
croscopic dynamical mechanism underlying this process 
has been first explored in [IT] , and has lead to the dis- 
cussion of charge asymmetry in large bodies [12] such 
as neutron stars, Galaxies and large scale cosmological 
magnetic field domains. 

Kharzeev et al [15] proposed to look for charge asym- 
metry induced by chiral CP-odd fluctuations at Rela- 
tivistic Heavy Ion Collider (RHIC), by means of charged 
pions' correlations. They pointed out that very intense 
magnetic field — up to 10 19 Gauss — is generated by the 
nuclei and also by motion of excited matter, which has 
certain overall positive charge density and vorticity, ex- 
cept in the case of exactly central collisions. The asym- 
metry in charged pion correlations relative to the collision 
plane has been recently observed by the Star Collabora- 
tion [13] [T3], see Fig.Q. The observed effect is indeed 
absent for central collisions. Still, the observed correla- 
tion is by itself CP-even ^< 9 2 > and in principle it may 
be induced by some other effect, absent in the usual event 
generators. To be sure that it is indeed proportional to 
< 9 2 > (averaging over the volume and events) some 
other experimental confirmation would be desirable. Di- 
viding it out would allow to single out the QGP response 
to 9, a property opening a window to new transport coef- 
ficients [El [17]. In this paper we propose CP-odd decays 
as such possibility. 




FIG. 1: A sketch of a chirally asymmetric spot (inner cir- 
cle) in a fireball at freezeout. The external magnetic field 
(not shown) is assumed to be directed vertically: induced 
electric field leads to charge separation (plusses upwards, mi- 
nuses downwards). The spot's location is biased upward, as 
explained in the text. During the decays, the ++ correlation 
is quenched by ambient matter less than H — . 



is in the hadronic phase, but the topological fluctuations 
have not washed out yet and there is a domain of finite 
size where the QCD vacuum is in a state with an effective 
9^0. In this regime one can use Chiral Perturbation 
Theory (ChPT), keeping in mind that the coefficients of 
the Chiral Lagrangian change their values with temper- 
ature . For our purposes we only need the CP-odd part 
of the ChPT Lagrangian in the leading order in 9. 

The ChPT Lagrangian is usually represented as a func- 

eit 11 ^] where n is 



tion of a non-linear field U (x) — 
the nonet of pseudoscalar mesons 
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The singlet and octet fields 0° and s are related to the 
physical rj, rj fields via the mixing angle olm as 



V 



ip s cos ctM — f° sin olm 
ip 8 sin ajf + y° cos olm- 



(3) 
(4) 



The CP-odd part of the ChPT Lagrangian in the leading 
chiral order and at leading order in the 1/N C expansion, 
is equal to [18] 

Leff = -*^f^~ [Tr[C/(x) - U\x)] - 21ogdet U(x)](5) 



II. CP FORBIDDEN DECAYS 

In this paper we are discussing the CP-violating pro- 
cesses in the phase of heavy ion collisions at the freezeout, 
which occurs at temperature T ~ 170MeV. The system 



where 9 is an effective 9 parameter. 

s FlN r ml „ 



(6) 



To study decay processes of a pseudoscalar meson into 
two pseudoscalar mesons, we then have to expand Eq. ^ 
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to third order in the meson fields and select all the inter- 
action terms which are compatible with a decay process, 
that is to say terms with a field of rest mass bigger than 
the sum of the rest masses of the two other fields. A sim- 
ple calculation reveals that the only decays allowed are 
77's into two pions 
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and we will only consider the process 77/77 — > 7r + 7r~ be- 
cause it is simpler to detect charged pions. The am- 
plitude of this process at leading order is given by the 
square of the coefficient of the corrispondent term in the 
Lagrangian 
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for value of the mixing angle we used the value olm — 
— 15°, see [19] . The topological susceptibility and the 
pion decay constant have to be substituted with their 
values at the freezout temperature Tj ~ 170MeV. 
The decay probability is then equal to 
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where in the last approximate equality we have substi- 
tuted the numerical values of the parameters involved, at 
T = 0. (We return below to the issue of T dependence). 

We now make a comparison to allowed decays: known 
decays into 77 — s- 37r are so much suppressed that the 
corresponding partial widths are about half of the total 



width of r„ = 1.18±.llfceV. So, with 
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needed to explain CP-odd fluctuations at RHIC, we actu- 
ally find that the forbidden decay rate in the asymmetric 
spots is orders of magnitude larger than the allowed one! 
Still, the corresponding lifetime is much larger than the 
duration of the freezeout Ary, after which all particles 



including etas leave the spot. Thus the fraction of etas 
to decay via the forbidden channel is still small 
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of the order of one per mill or so. (The rest decay in the 
usual way outside the spot.) 

We now discuss how masses and parameters of the chi- 
ral Lagrangian (like topological susceptibiltics) are modi- 
fied from their vacuum values, at finite temperature. We 
start by considering the mass of eta, which is given by 
Gell-Mann-Oakes-Renner relation and thus 



< ss> (T) 

1WT' 



(13) 



While the quark condesate changes appriciably for T — > 
Tj , the same is true for /„■ , f v , and the ratio is believed 
to be changed by very small amount. We looked at avail- 
able lattice literature, and although there are measure- 
ments of all three quantities m^(T),< ss > (T),f„(T) 
we have not found common datasets in which all can be 
compared together. Lacking that, uncertainties due to 
different lattices and quark masses are too large to give 
some particular numbers. Yet nothing prevents to do so 
in a dedicated calculation. 

The topological susceptibility does not change appre- 
ciabily and account just for a small factor in the eta mass: 
it may however change T,^^ and - more importantly - 
the value of the 77' mass. 

We will now subsequently discuss the factors which af- 
fect the invariant mass window in which decays are to be 
observed. These include: 

(i) spatial dependence of effective 9 parameter 

(ii) meson mass modification as a function of tempera- 
ture. 

(i) The spatial dependence of the effective parame- 
ter 9{x) can be viewed as additional 3-momentum k w 
1 / R sp ot coming from its Fourier transform. (We do not 
include the time-dependent component, assuming that 
diffusion of chirality is a slow process, and the corre- 
spondig frequency is negligible.) As a result, the invari- 
ant mass of the final state is modified by 



m 2 — k 2 



2m 2 



(14) 



The magnitude of the spot size and the corresponding 
momenta will be discussed in the next section, for es- 
timate we will now use k — 50 — 100 MeT^. From the 
formula above we see that the correction to the mass is 
then in the range of -(0.005 ~ 0.02). 

(ii) thermal mass modification of the eta, can also be 
viewed as the energy shift of a particles sitting in a cer- 
tain potential created by all other memebers of the fire- 
ball, at the moment of the decay. In our case we cosider 
two different inital particles, 77, 77', which are related with 
different phenomena and thus should experience different 
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shifts. The former one is mostly Goldstone meson (ignor- 
ing their mixing), the SU(3) partner of the pions, thus 
its mass is obtained from Gell-Mann-Oaks-Renner rela- 
tion; the latter, 77 , has a mass given by Witten-Veneziano 
formula. While the quark condensate and the decay con- 
stant are expected to have significant changes themselves, 
and in the chiral limit vanish at T = T c , the ratio appear- 
ing in the eta mass is believed to change with T much less. 
For example, the quark condensate at T c h ~ 170 MeV is 
reduced by about 20% [20]. 

At the end of this section let us briefly discuss the re- 
laxation of the spot itself. During the freezeout period, 
when the CP-odd spot is populated with mesons, some 
energy/momentum exchange between 9(t, x) and these 
mesons is possible and is described by the same chiral 
Lagrangian: for example, it was included in our discus- 
sion of the visible width of the forbidden decays above. 
What happens after that period, when the residual 9(t, x) 
remains alone, as a vacuum perturbation? Depending on 
the energy available, it can still decay into a number of 
channels. In the order of reduced mass those would be 
e.g. processes like 9 — > 77', or 9 — > i]tttt or eventually 
9 — > 77 till it relaxes completely to 9 = 0. The rate of 
all of them, if needed, can be calculated from the same 
standard chiral Lagrangian. 
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2: Preliminary STAR data on charge correlation [131 



III. WHAT IS THE SIZE OF THE "CHIRALLY 
ASYMMETRIC SPOT"? 

Let us return to the geometry of the charge correlations 
shown in fig|2] Note that we placed the asymmetric spot 
at the periphery of the fireball: we did so for the following 
reason (known as trigger bias): such geometry increases 
chances that both charge pions (say 2-7r + ) are observed, 
escaping from the system unscattered and moving in a 
certain direction (up in the figure). There is no such bias 
effect for opposite charges H — since those are expected 
to move approximately in opposite direction. Existing 
experience with jets show that one of the two or both 
would be rescattering and thus quenched by matter much 
more. So, the observed ratio of opposite-to-same charge 
effects R(-\ — / + +) is expected to be small. 

Now look at the STAR data shown in figj2) Those 
for AuAu collisions (shown by the closed points) indeed 
finds the effect for opposite charge -I — correlation to be 

much smaller than for same charge + + / one, R(-\ — 

/ + +) -C 1. However more recent data taken for CuCu 
collisions (open points) show a nearly symmetric picture, 
with \R(-\ — / + +)| « 1 (up to a sign, as expected). We 
take it as the direct experimental evidence for the fireball 
sizes at freezeout having the following relation to the size 
of the CP-odd spot 

RAuAu ^ Rspot ~ RcuCu (15) 

We used this last approximate equality for the estimate 
of the corresponding momentum in the previous secton. 



IV. SUMMARY AND DISCUSSION 

The charge pion correlations observed by STAR col- 
laboraton [13] [H] is by itself CP-even. To be sure that 
it is indeed proportional to < 9 2 >, we propose to look 
for CP-forbidden decays such as 77 — > 7r + 7r~, 77' — * 7r + 7r~. 
According to our estimates, the rates of the decays at 
freezeout are orders of magnitude larger than the allowed 
3-meson modes, and the total fraction of the forbidden 
decays should be of the order of 10~ 3 or so. 

We think that with the available statistics of ~ 10 9 
events one can separate the peaks in the invariant mass 
from the statistical background of random Tt + ir~ pairs 
with the same invariant mass. (Since the channels are 
very different, the usual decays make no background.) 
From theoretical point of view 77 stands a better chance 
because its mass is expected to change little, as it is pro- 
portional to the square root of the ratio of quark conden- 
sate and F w , with similar T-dcpcndcnce. The 7/ mass 
is related to topological susceptibility and may change 
more. But from experimental point of view, due to the 
values of 77 and 77 rest mass, the signal/background ratio 
may be larger for 77 decays: since the freezeout tempera- 
ture it is Tf ~ YlQMeV , RHIC collisions are more likely 
to produce pions with momentum of the order of 0.5 
than 0.5 m^. The issue of course require further studies, 
including experimental resolution and acceptance: but, 
we repeat, it does not require new data, just their anal- 
ysis. 

If observed, the CP-forbidden decays would provide 
an independent measure of the CP-odd fluctuations. Di- 
viding < 9 2 > out would allow to single out the QGP 
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